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Long-range charge transfer excited states are notoriously badly underestimated in time-dependent 
density functional theory (TDDFT) . We resolve how exact TDDFT captures charge transfer between 
open-shell species: in particular the role of the step in the ground-state potential, and the severe 
frequency-dependence of the exchange-correlation kernel. An expression for the latter is derived, 
that becomes exact in the limit that the charge-transfer excitations are well-separated from other 
excitations. The exchange-correlation kernel has the task of undoing the static correlation in the 
ground state introduced by the step, in order to accurately recover the physical charge-transfer 
states. 
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I. INTRODUCTION 

How or whether excited states of long-range charge 
transfer character are captured in time-dependent den- 
sity functional theory (TDDFT) has recently received 
much attention [1-9], partly because of its importance 
in systems of biological and chemical interest. Many 
of these molecules are big enough that traditional wave- 
function methods to calculate excited states become pro- 
hibitively expensive; TDDFT, on the other hand, scales 
favorably with the number of electrons, while remaining 
reliably accurate for many excitations. There is therefore 
much interest in applying TDDFT to such systems, and 
understanding why, when, and where it is expected to 
work well. 

TDDFT is an exact theory, based on the mapping 
of the interacting electronic system to a non-interacting 
one [10], far faster to solve. In practise the unknown 
time-dependent exchange-correlation potential must be 
approximated, introducing errors in excitation energies. 
These are not always well-understood, although there 
has been some recent progress [11], for example in the 
significance of the asymptotic behaviour of the poten- 
tial [12, 13], and in the severe frequency-dependence 
needed for double-excitations [14-16]). 

It has been found [1-9], that excitations involv- 
ing charge-transfer (CT) between widely separated 
species within a molecule are severely underestimated in 
TDDFT. This is in contrast to most excitations, which 
are accurate to within a few tenths of an electron-volt. 
This has serious consequences. For example [3], predic- 
tion of CT quenching of fluorescence in light-harvesting 
bacteria, contrary to observation. 

There have been several recent attempts to overcome 
this charge-transfer problem. In Ref. [8], a ASCF cor- 
rection is made, while in Ref. [1], it is proposed that 
TDDFT be mixed with configuration interaction singles. 



'Electronic address: nmaitraOhunter.cuny.edu 



In Ref. [9] , a long-range correction using Hartree-Fock is 
implemented. Ref. [2] suggests a simple improvement by 
shifting the donor HOMO and acceptor LUMO energies 
based on the derivative discontinuity [17]. Ref. [4], on 
the other hand, remains purely within TDDFT, and sug- 
gests an empirical asymptotically-corrected kernel that 
diverges with separation. 

In the present paper, we also stay strictly within 
TDDFT, and ask, in contrast to all other approaches, 
what is the exact TDDFT description of charge-transfer 
in one simple case: charge-transfer between two open 
shell species, when the charge transfer states are well- 
separated from the other excitations of the system. By 
solving a simple model we deduce what the exchange- 
correlation kernel must be in this case. 

An important feature of the present analysis involves 
the step in the exact ground-state Kohn-Sham (KS) po- 
tential that appears between two widely-separated open 
shell species of different ionization potentials. This step 
is also present in orbital approximations such as exact 
exchange (or the KLI approximation to it [18]). It is not 
present in local or semi-local gradient approximations: 
so our model for the exact exchange-correlation kernel 
does not apply to these cases. Our purpose here is to ex- 
plore how charge-transfer is described exactly in a simple 
model case. 

We begin in Section II with a review of TDDFT lin- 
ear response theory, and present the problem of long- 
range charge transfer. A simple model to study this is 
introduced in section III. We describe the significance 
of the step that develops in the KS potential when a 
(closed-shell) molecule composed of two different open- 
shell species are pulled apart. The KS ground-state is of 
a fundamentally different nature than the true ground- 
state, resembling the problems of homonuclear disso- 
ciation; the KS charge-transfer energies become zero! 
We show how inclusion of the electron-electron interac- 
tion breaks the degeneracy and yields the correct nature 
of the ground-state and the CT states. In section IV 
we describe the implication of static correlation for the 
TDDFT exchange-correlation kernel, and derive a model 
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for it that is valid in the limit that the charge-transfer 
excitations are well-separated from all other excitations 
in the system. The features are strong-frequency depen- 
dence and exponential dependence on the separation of 
the two species. Finally, in section V, we discuss other 
ways that static correlation haunts TDDFT, so far unex- 
plored. 



II. TDDFT LINEAR RESPONSE AND CT 
STATES 

We begin by briefly reviewing the TDDFT linear re- 
sponse formalism. 

Although the density of the KS system is defined to 
be that of the interacting system, the excitation ener- 
gies are not the same: linear response theory tells us 
how to correct them. Applying a small perturbing po- 
tential to a ground-state, and measuring the density re- 
sponse defines the susceptibility or density-density re- 
sponse function, x[no](r,r',t-t') = Sv ext (rt)/Sn(r't')\ no . 
The susceptibility of the true interacting system is related 
to that of its non-interacting Kohn-Sham counterpart, 
Xs[«o](r,r',t-i') = 6v s (rt)/5n(r't')\ no , through an inte- 
gral equation, written in the frequency-domain as [19]: 



(i) 



Here the Hartree-exchange-correlation kernel is the sum 
/ H xc[™o](r,r',w) = / H (r,r') + / xc [n ](r, r', w). The 
Hartree kernel is the density-functional-derivative of the 
Hartree potential, / H (r,r') = l/|r— r'|, and the exchange- 
correlation kernel is that of the exchange correlation po- 
tential, f xc [no](r, r', t — t') = 5v xc (rt) / 5n(r' , t'). Tran- 
sition frequencies of the true system lie at the poles of 
x(r, r', u>), and oscillator strengths may be obtained from 
the residues. The poles of x s (r, r', w) are at the KS sin- 
gle excitations; these are shifted to the true excitations 
through the action of the Hartree-exchange-correlation 
kernel. So, Eq. (1) enables us to obtain the interacting 
excitation energies and oscillator strengths from the KS 
susceptibility and the Hartree-exchange-correlation ker- 
nel. In principle, the exact spectrum of the interacting 
system is obtained; in practise, approximations must be 
made for (a) the xc contribution to the ground-state KS 
potential, and (b) the xc kernel / X c(^)- 
In the Lehman representation, 



X(r,r';w) = ^ 



F f (r)i?(r') f;(r)F 7 (rO ] 
+ i0+ J 



(2) 



where, Fj(r) — (0|n(r)|7) with n(r) being the one-body 
density operator, / labels the excited states of the inter- 
acting system, and u>i is their transition frequency. This 
expression also holds for the KS susceptibility where the 
excited-states are excited Slater determinants and the 
transition frequencies are orbital energy differences. Ex- 
citations of atoms and molecules are often calculated in 
a matrix formulation of these equations [20]: one lets 



q = (i,a) be an index representing a single excitation: 
a transition from an occupied KS orbital (f>i to an un- 
occupied one <p a , and let ui q be the difference in the KS 



orbital energies, uj q 



6i. Then, the squares of the 



true transition frequencies f2/ 
of the matrix 



Loj are the eigenvalues 



Sl(u) qq > = 5 qq >uJ 2 q + 4i/Lj q Lu q >[q\f HXC (u)\q'] 



(3) 



where 



[91/hxcW] = J drdrV*(r)0a(r)/ H xc(r,r'^)^(r')^(r'). 

(4) 

Oscillator strengths of the true system are related to the 
eigenvectors [20]. 

When the coupling between excitations is very small, 
one may neglect the off-diagonal elements of the matrix, 
producing the "small-matrix approximation" [11, 21, 22]: 



tu 2 q + 4w 9 [g|/Hxc(w g ) 



(5) 



This is equivalent to keeping only the backward and for- 
ward transitions at frequency ui q in the Lehman repre- 
sentation of the response function. When the shift from 
the KS transition is small, one may simplify this further 
to get the "single-pole approximation" [19]: 



uj = u q + 2[q\f KXO (uj q )\q] . 



(6) 



Now, consider applying the theory to long-range 
charge-transfer excited states. In the limit of large sepa- 
rations R, the exact energy cost for transferring an elec- 
tron from donor to acceptor, is 



ui = I D - A A - l/R 



(7) 



where Id is the ionization energy of the donor, and Aa 
is the electron affinity of the acceptor. As discussed in 
Refs. [1, 2], the failure of TDDFT to reproduce this is 
evident from a single-pole analysis (Eq. (6)), where 



uj = e A - e 



g + J d 3 rid 3 r 2 F(ri)/ H xo(r ) r' ) a;,)J , (r 2 ) 
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= I D -A S 



(8) 



Here F(r) = 0g(r)<^(r) is thc product of the HOMO 
of the donor (KS orbital energy e^) and the LUMO of 
the acceptor (KS orbital energy e^). Because there is 
exponentially small overlap between the atomic orbitals 
on the widely separated donor and acceptor, the integral 
term vanishes, and the TDDFT energy collapses to the 
bare KS energy, as indicated in the second line. We have 
used the fact that the KS HOMO energy is exactly the 
negative of the ionization energy, but the LUMO energy 
differs from the negative affinity by the discontinuity [17, 
23-25]: 



A 



Ay 



-e L + A* 



(9) 
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Now, common approximations underestimate /, but how 
does exact TDDFT get the exact energy? That is, sup- 
pose we had the exact ground-state Kohn-Sham potential 
and also the exact exchange-correlation kernel. Then the 
exact TDDFT HOMO energy is indeed I D , but how docs 
the exact TDDFT retrieve the discontinuity A xc , and the 
— 1/R Coulomb fall-off ? We address this in the next two 
sections. 



III. CT BETWEEN TWO OPEN SHELL 
SPECIES: THE ROLE OF STATIC 
CORRELATION 

To study this question, we consider first the simplest 
model: electron transfer between two onc-clcctron neu- 
tral "atoms", separated by a large distance R. Let ii(2) 
be the ionization energy of atom 1(2); then the atomic 
orbital occupied in the ground-state has energy e H = 
—I\{2), where the H subscript stands for HOMO of the 
atom. When we consider the closed-shell molecule com- 
posed of the two atoms at large separation, the ground- 
state KS potential develops a step between the atoms, 
that exactly lines up the atomic HOMOs [23, 26]. To see 
this, consider the Kohn-Sham wavefunction $o for the 
molecule. This is the doubly-occupied bonding orbital: 

*o - 0o(ri)0o(r 2 )(|n>-|IT»/V2, with 

-Mr) = (0i(r) + fc(r))/V2(l + Si 2 ) (10) 

where Si 2 is the overlap integral J 0i(r)</> 2 (r)c? 3 r, expo- 
nentially small in the separation R. In the limit of large 
separation, 0i j2 denotes the occupied atomic orbital of 
atom 1(2). Inserting tfi into the KS equation, while rec- 
ognizing that near each atom, the orbital must reduce to 
that atom's atomic HOMO, with the appropriate orbital 
energy, one sees that the exchange-correlation potential 
must develop a step of size \I 2 — h\ between the atoms; 
raising the potential of the atom with the larger ioniza- 
tion energy. A simple 1-d example is given in Fig. 1 where 
each "atomic nucleus" is an Eckart potential. The poten- 
tial eventually goes back down to zero on the right-hand- 
side. (A similar picture for two delta-function "nuclei" 
can be found in Ref. [23] ) . The step renders the HOMOs 
of each atom degenerate, and is necessary to prevent dis- 
sociation into fractionally charged species [17]. If they 
were not exactly lined up, then one could lower the en- 
ergy of the system by transferring a fraction of charge 
from one of the atoms to the other, contradicting the 
charge neutrality that is observed in nature. This phe- 
nomenon is closely related to the derivative discontinuity 
discovered in the 80's [17, 23-26]. It is important to 
note that the step is present in exact DFT and is ap- 
proximated in spatially-nonlocal approximations such as 
exact-exchange [18, 27] Most common approximate po- 
tentials, such as LDA, garner only local or semi-local 
density information, and so do not display the step. 

The fact that in exact DFT the atomic HOMO's be- 
come degenerate has an apparently drastic consequence 
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FIG. 1: The 2-electron KS potential for two widely separated 
"Eckart" atoms, Kxt(a') = l/cosh 2 (a;) + 1.8/cosh 2 (a; - 15). 
The highest occupied orbitals of each atom are shown solid, 
and the first excited state of the right-hand atom is shown in 
dash-dot. 

for CT excited states: Whereas a (semi-)local DFT ap- 
proximation would yield a finite KS transition energy, 
albeit underestimated, the exact KS system yields zero 
energy for charge transfer between the atoms! 

The Kohn-Sham energies are poor because, in this 
case, the KS description is fundamentally different from 
the interacting system. The KS ground-state involves 
each electron delocalized over both atoms, whereas the 
true wavefunction has one atom on each. We shall come 
back to this point shortly. 

The alignment of the atomic HOMO's leads to three 
nearly degenerate singlet two-electron KS states, whose 
spatial parts are: 

$o = 0o(ri)<A)(r2) 

$ 9 = (0 o (ri)0o(r 2 ) +0o(ri)0 o (r 2 )) /V2 
$ D = 0o(ri)0 o (r 2 ) (11) 

where 4>o is the bonding orbital of Eq. (10), and <po is the 
antibonding orbital 

0o= (&-02)A/2(i -S 12 ). (12) 

At large separations the "singly-excited" $ g is slightly 
higher in energy than <f>o by tunnel splittings between 
the atomic orbitals, and the "doubly-excited" <J>d higher 
again; the orbital energy of the antibonding state is at 
w ~ exp(— R), where R is the separation between the 
atoms. The proximity of $q and <&d to $o is a signature 
of the static correlation in the ground-state. 

Before discussing the implication of this static 
correlation for TDDFT, wc first show that the elec- 
tron interaction splits the degeneracy between the 
Kohn-Sham states, recovering states of CT nature and 
approximations for their energies. We assume that the 
tunnel-splittings are much smaller than the energy- 
splitting from the electron-electron interaction, and that 
the basis of the three states, Eqs. (11), has negligible 
couplings to all other KS excitations. To simplify the 
discussion, we choose / 2 > I\. Diagonalizing the true 
Hamiltonian in this basis, yields the three states: 
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(i) the Heitler-London ground-state = 
(<Mr)0 2 (r') + 2 (r)<Mr'))V2 

(ii) charge-transfer excited state from atom 2 to atom 
1: \?2-»i = (j>i (r)0i (r'). This has energy, relative to the 
Heitler-London state, 

u) X =h-A 1 - l/R (13) 

where A = A s + A^' ox , with 

A~ = - J d 3 r y"dV^(r) 2 ff (r') 2 K=o(r-r') (14) 

Here V ee is the electron-electron interaction, V ee = l/l r — 
r'|. For later purposes, the subscript H stands for the 
highest occupied orbital of the atom; for one-electron 
atoms, it simply indicates the occupied atomic orbital 
on the atom. Note also that for open-shell atoms, A s — 
I=-e H . 

(iii) charge-transfer excited state from atom 1 to atom 2: 
*i^2 = 02(r)02(r'), of energy 

uj 2 = h-A 2 - l/R (15) 

Thus, the correct nature of the ground-state, and CT 
states with the correct form of energy (c.f. Eq. 7), in- 
cluding an approximation for the exchange-correlation 
contribution to the electron affinity, are recovered. 

In a sense, accounting for the electron-interaction un- 
does the static correlation in the system. Because of the 
step in the potential, and the consequent nearly degener- 
ate triple of KS determinants, the ground-state problem 
resembles the well-known homonuclear dissociation prob- 
lems of ground-state DFT [28, 29], where the dissociation 
limit is incorrect because the KS ground-state involves 
mixing ionic combinations into the true neutral Heitler- 
London form. Static correlation was an important agent 
above in treating the perturbation: The Hamiltonian can 
always be written in terms of the KS Hamiltonian H s , as 
H = H s + V co — v H — v X c- For any two-electron sys- 
tem, as we have here, w x = —Vn/2 ~ O(A), where A is 
the interaction strength. Typically the correlation poten- 
tial is 0(A 2 ), appearing at the next order of interaction 
strength, but when there is static correlation, it is much 
stronger: v c cancels v H + w x and adds the step, yielding 

D HX c( r ) = , for r near atom 1 

= I2 — I\ , for r near atom 2 

This must be so in order for the KS equation for the 
molecular orbital 0o to reduce to that for the atomic 
orbital 01,2, respectively, near atom 1 and atom 2. 
We discuss how static correlation affects the exchange- 
correlation kernel of TDDFT in the next section. 

We now consider the approximate result for A xc , 
Eq. (14). The derivation neglected all KS states ex- 
cept for the three in the basis (11): the approximation 
thus becomes exact when the basis is truly isolated, or in 
the limit of weak interaction so that the higher-lying KS 
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FIG. 2: Exchange-correlation contribution to the electron 
affinity of a fermion in a 5-well: exact (extracted from 
Ref. [30]) and approximate (Eq. 14) 

states are not appreciably mixed in. A simple example 
for which the exact affinity at all interaction strengths is 
easily calculable, is that of a fermion in one dimension 
living in a delta-function well (i.e. a one-dimensional H 
atom). It is a standard textbook problem to find the 
only bound-state of the well. The electron affinity can 
be obtained by subtracting its ground-state energy from 
that of two interacting fermions in the well. When their 
interaction is a delta-function repulsion of strength A: 

H=-jz\^2-jZ 5 ^) + X8 ^- x ^^ ( 16 ) 

i 1 i 

the ground-state energy has an exact solution [30]. 
(Atomic units are used throughout the paper). In fig- 
ure 2, the exact A xc has been plotted from ^4 XC = 
E(l) - E{2) -A s = -E{2) - 21, where the ground-state 
energy for two fermions, E(2), is obtained from Ref. [30], 
and that for one fermion is £'(1) = — ^4 S = ~I = —0.5H. 
Comparing this with the approximation from Eq. (14), 
where 4>h is the textbook solution for the orbital in the 
one-fermion problem, we see that, as expected, the ap- 
proximation becomes exact in the weak interaction limit. 

This result for A xc is consistent with a perturbative 
analysis of E(N) — E(N + 1). Consider adding one elec- 
tron to an N — 1-species. In an unrelaxed approxima- 
tion, the extra electron plops into the same spatial orbital 
as the first. A perturbative approximation for the elec- 
tron affinity is obtained from calculating the difference 
in the Hamiltonians of the 2-electron and 1-electron sys- 
tem, in the unrelaxed doubly-occupied orbital. One finds 
exactly the result, Eq. (14). 

Although our results so far have been derived for 
charge transfer between two one-electron species, it 
is straightforward to show that they hold for charge- 
transfer between two general open-shell species. There 
the KS molecular HOMO sits atop a spin-paired KS de- 
terminant, containing all the inner orbitals of both atoms. 

So far, we have examined the exact ground-state KS 
potential for a closed shell molecule composed of two 
open-shell species, and found that the step in between the 
species leads to a near zero value for charge-transfer in 
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the Kohn-Sham system. We showed that the static corre- 
lation is broken by the electron-electron interaction, and 
how the latter allows us to recapture the true ground- 
state and excited charge-transfer states, together with 
an approximation for the exchange-correlation electron 
affinity. 

IV. THE EXCHANGE-CORRELATION 
KERNEL FOR CHARGE-TRANSFER 

How does the static correlation in the ground-state 
affect the TDDFT description of charge-transfer? In 
TDDFT, the exchange-correlation kernel must play the 
role of the diagonalization of the previous section in sal- 
vaging the CT energies from their near zero KS value. 
We now examine what the nature of the exact f KC must 
be, and the signatures of static correlation in the kernel. 

Consider the KS response function, Xs( r , r',u). We 
may write this as 

Xs (r,r»« 2 2 " X s (r,r',(u)) (17) 

to — [CO ) 

where we include only the contribution to Eq. (2) 
from the antibonding transition, valid in the limit 
that the KS bonding-antibonding pair is well-separated 
from the higher KS transitions. The assumption 
is that to lowest order in the interaction strength, 
the interacting charge-transfer states arise from this 
transition only. This is analogous to the ap- 
proximation of restricting to the 3x3 subspace 
when we considered diagonalization. The numerator 
X s (r,r', H) = (4>?(r) - <%(r)) (<t>\(r>) - ^(r')) /4 + 
O ((uj 2 - Q 2 )A s /[Q(lo 2 - Ag)]) , where A s is the fre- 
quency of the next highest KS transition. 

Notice that %s vanishes exponentially with separation 
of the atoms due to the exponential dependence of the 
tunnel splitting uj on separation. This can be understood 
from x s = j^(lu): the perturbation at Co excites the 
antibonding transition, whose density difference from the 
bonding orbital's is exponentially small with separation. 

Now we turn to the interacting response function, 
X(r,r',w), which is a sum over all excitations of the in- 
teracting system. As in the KS case, we zoom into the 
excitations that describe charge transfer: 

*( r > r ' w ) ~ 772 — ^2 X i( r > r > H)+7T5 — 772 x 2(r, r , (w)) 

UJ — UJ-^ UJ — 

(18) 

Within the isolated basis approximation, Xi(r, r', (oj)) = 
X 2 (r,r',(w)) = 20i(r)0 2 (r)0i(r / )0 2 (r / ) + 0(1/ A). 
Again the response function vanishes exponentially with 
separation, but this time because the CT transitions have 
exponentially weak oscillator strength due to overlap be- 
tween widely separated spatial regions: in the orbital 
product 0i(r)0 2 (r), each orbital is exponentially small as 
a function of atomic separation R in the places where the 
other is finite. The inverse functions, x an d Xs> are 



therefore both exponentially large as a function of sepa- 
ration R (not to be confused with the finite behaviour as 
a function of |r — r'|). 

We note that we are retaining backward transitions in 
the expressions above: these are "small-matrix approxi- 
mations" rather than "singlc-polc approximations" that 
include only the forward transition. The latter is not 
valid here because the antibonding transition frequency is 
very small: the backward transition is of almost the same 
magnitude as the forward, so should not be neglected. 
The single-pole approximation is only valid when the 
shift from the KS energy is small [11]: here, the shift 
is the entire CT energy, since the KS transition energy is 
near zero (i.e. exponentially small in the separation). 

As we saw in the diagonalization procedure, the CT 
transitions of the interacting system arise from the KS 
subspace of bonding and antibonding orbitals. In the 
diagonalization procedure, the number of states was pre- 
served: the three KS 2-electron singlet molecular deter- 
minants of Eq. (11), composed of the bonding and anti- 
bonding orbitals, were rotated by the electron-electron 
interaction into the Heitler-London ground-state and CT 
states onto each of the two atoms. The picture is some- 
what different in the response functions. Here the space 
is of excitations out of the ground-state, so in the inter- 
acting case there are two (one to each CT state). But 
the KS response function has only the single-excitation 
(corresponding to <I>o — » ® q of Eq. (11)): double excita- 
tions cannot appear in the KS response function because 
the numerator involves the one-body density operator, 
which connects only states differing in one orbital (see 
also Ref. [14]). It is the job of the exchange-correlation 
kernel, / xc to generate an extra pole, and also to mix 
them. Specifically, 

/hxoM = x s T 1 M-x~ 1 M 

_ 1 f ^-oj 2 ) (u^-u^-u^ A 

2\ u s uj 2 (uj 2 - uj 2 ) + uj^oj 2 - uo 2 ) 1 J 

Due to the /-sum rule, in the limit of the isolated basis, 

loX s = wiXi + uj 2 X 2 « (wi + uj 2 )X 1 (20) 

we find X7 1 = (u)\ + ui-zjX^ 1 /u>. Finally, we require 
that a dressed small-matrix approximation, arising from 
plugging Eq. (19) into Eq. (5) yields the true transition 
frequencies u> = us\ and u> = ui 2 as solutions. This fixes 
[glAjT 1 ^] = 1/2, and our expression for the kernel be- 
comes: 

^M/ HXC M| ? ] = .L + ^^ + 4^ (2D 

4 U) — U1\L0 2 

for the KS transition q =bonding^anti-bonding orbital, 
and where ui av = (u>i +lo 2 )/2. This, with u>i and ui 2 given 
by Eqs. (13), (14), and (15), gives the exact exchange- 
correlation kernel matrix element in terms of KS quanti- 
ties, in the limit that the charge-transfer states are iso- 
lated from the other transitions, in that coupling to them 
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can be neglected. Notice the strong non-adiabaticity in 
the exact kernel, manifest by the pole in the denominator 
on the right-hand-side. 

In the usual presentations of the TDDFT charge- 
transfer problem, the failure is due to the vanishing over- 
lap between the occupied orbital, on one atom, and the 
unoccupied one on the other atom, widely separated. 
This means that q = <fii(f> a vanishes exponentially with 
separation. This does not occur in the exact analysis 
of CT between open shells presented here: here the oc- 
cupied orbital is the bonding orbital and the unoccu- 
pied is the anti-bonding orbital, so their overlap goes as 
(4>i - <t>l) / 2 where fa is the atom l's HOMO (or LUMO) 
and fa is atom 2's HOMO (or LUMO). Although the 
overlap remains finite, the matrix element [<ji/ H xc|<z] di- 
verges exponentially with interatomic separation R, as 
1/u), as can be seen from Eq. (21). 

A divergence of / HX c-matrix elements has previously 
been found in H 2 dissociation [29], where the the lowest 
singlet-singlet KS transition energy vanishes with sepa- 
ration, while the true energy approaches a finite num- 
ber; the step in the KS potential for our hcteroatomic 
molecules leads to a similar feature raising its head again 
here, but now with the kernel being strongly frequency- 
dependent. Also, the CT treatment in Rcf. [4] utilized an 
empirically determined kernel that displays also exponen- 
tially large behaviour as a function of atomic separation. 

V. DISCUSSION AND OTHER IMPLICATIONS 
OF STATIC CORRELATION 

By considering electron transfer between two one- 
electron species at long-range, we investigated how 
charge-transfer is captured in the exact KS system. The 
bare KS energies for CT approach zero exponentially 
with the separation of the species, but including the 
electron-electron interaction, splits the near-degeneracy 
and recaptures finite CT energies of the correct form. 
We derived the form of the TDDFT exchange-correlation 
kernel, Eq. (21), that becomes exact in the limit that the 



charge-transfer states are isolated from all the other ex- 
citations in the system. The main features of the kernel 
are (i) a strong frequency-dependence, due to the mixing 
in of the double excitation to the antibonding KS state, 
and (ii) a dependence on the inverse tunnel splitting be- 
tween the two species, that goes exponentially with their 
separation. 

A crucial feature of the exact KS system in this paper is 
the step between the two widely separated species. This 
gives rise to static correlation in the ground-state, and 
hence the strong frequency-dependence and exponential- 
dependence on the atomic separation of the / xc matrix el- 
ement, discussed above. Not only does it play a vital role 
in the description of CT states, but it also has very inter- 
esting, and unexplored, effects on other excitations of a 
long-range heteroatomic molecule. For example, higher 
atomic excitations become KS resonances, as illustrated 
in Figure 1. Here two one-dimensional "Eckart atoms" 
are shown: the second excited state of the atom on the 
right can tunnel out of the barrier presented by the step. 
The acrobatics the exchange-correlation kernel must per- 
form in order to turn such a resonance back into a bound 
state of the true system, will be pursued in future work. 

Another interesting consequence of the step is that for 
every single excitation out of the HOMO bonding orbital, 
fa — > 4>a, there is a nearly degenerate double-excitation, 
(fa, fa) — > (fa fa), where the other electron occupying 
the HOMO is excited to the antibonding orbital, sepa- 
rated in energy only by the tunnel-splitting. Doubles are 
absent in the KS response [14], but they are essential in 
this case for a correct prediction of the nature and energy 
of the interacting state: for example, it is needed to lead 
to a charge-neutral excited state on each atom [31]. The 
presence of double excitations immediately leads to poles 
in the frequency-dependence of the exchange-correlation 
kernel [14]; the static correlation implies the poles are 
ubiquitous and will be investigated in future work. 
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